課題４成果物デモンストレーション
Contents
#!pip install qiskit
#!pip install matplotlib
#!pip install pylatexenc
import qiskit
from qiskit import assemble, transpile
from qiskit import QuantumCircuit, ClassicalRegister, QuantumRegister
import qiskit.ignis.verification.quantum_volume as qv
---------------------------------------------------------------------------
ModuleNotFoundError Traceback (most recent call last)
Input In [2], in <cell line: 1>()
----> 1 import qiskit
2 from qiskit import assemble, transpile
3 from qiskit import QuantumCircuit, ClassicalRegister, QuantumRegister
ModuleNotFoundError: No module named 'qiskit'
# qubit_lists: QV回路を生成するための量子ビットサブセットのリストのリスト
nqubits = 10
qubit_lists = [[i for i in range(nqubits)]]
# ntrials: サブセットごとに作成するランダム回路の数
ntrials = 1
import warnings
warnings.filterwarnings('ignore')
qv_circs, qv_circs_nomeas = qv.qv_circuits(qubit_lists, ntrials)
qv_circs_nomeas[0] = qiskit.compiler.transpile(qv_circs_nomeas[0], basis_gates=['u1','u2','u3','cx'])
qv_circs_nomeas[0][0].draw(output='mpl')
f = open('hoge.qasm', 'w')
f.write(qv_circs_nomeas[0][0].qasm())
f.close()
課題４成果物デモンストレーション¶
・量子コンピュータ、初心者のための教育用GUI https://qulacs-gui.github.io/qulacs-simulator/¶
・高速シミュレータQulacsを用いたシミュレータバックエンド¶
・量子コンピュータ実機制約を満たしたQASMへと回路変換するトランスパイラ¶
Qulacsを用いたシミュレータバックエンド¶
シミュレータバックエンドの種類
ideal simulator with sampling
ideal calculation of observable
estimation of observable with ideal sampling (with shot noise)
noise simulator (sampling, observable calculation, observable estimation with sampling)
f = open('hoge.qasm', 'r')
hoge_qasm = f.read()
f.close()
print(hoge_qasm)
OPENQASM 2.0;
include "qelib1.inc";
qreg q[10];
u3(1.5469288,-2.8277906,2.7270131) q[0];
u3(2.2813222,1.4816381,-1.7218133) q[1];
u3(1.3242354,-1.7983311,0.541192) q[2];
u3(1.6351298,-2.2318667,3.1028801) q[3];
cx q[0],q[3];
u3(0.46365387,-pi/2,pi/2) q[0];
u3(1.5699358,-3.07503,-3.1286844) q[3];
cx q[0],q[3];
u3(0.03025315,0,-pi/2) q[0];
u3(1.5043758,0.0044345716,-1.6375118) q[3];
cx q[0],q[3];
u3(1.7285025,-0.19931365,-0.55479288) q[0];
u3(1.5249933,0.91057055,-3.0124968) q[3];
u3(0.44889372,-2.411138,1.4739187) q[4];
u3(1.5219969,0.5122379,0.48454526) q[5];
cx q[5],q[1];
u3(1.5929573,-3.0788164,2.8026727) q[1];
u3(0.72010473,-pi/2,pi/2) q[5];
cx q[5],q[1];
u3(1.5043758,0.0044345716,-1.6375118) q[1];
u3(0.089450441,0,-pi/2) q[5];
cx q[5],q[1];
u3(1.2805837,-1.0213828,0.69365036) q[1];
u3(0.80533913,1.7392758,-2.9558717) q[5];
cx q[3],q[5];
u3(0.59726811,-pi/2,pi/2) q[3];
u3(1.5790927,-3.0755427,3.0168207) q[5];
cx q[3],q[5];
u3(0.22036353,0,-pi/2) q[3];
u3(1.5043758,0.0044345716,-1.6375118) q[5];
cx q[3],q[5];
u3(1.5324398,1.4495826,-1.6591527) q[3];
u3(1.3084482,-1.5742808,-0.35063457) q[5];
u3(1.4359933,-1.882699,1.6719836) q[6];
u3(1.976004,2.6982712,-2.2666844) q[7];
cx q[7],q[4];
u3(1.6177998,-3.0944373,2.3583626) q[4];
u3(0.57344337,-pi/2,pi/2) q[7];
cx q[7],q[4];
u3(1.5043758,0.0044345716,-1.6375118) q[4];
u3(0.28823622,0,-pi/2) q[7];
cx q[7],q[4];
u3(0.92801429,-0.24644449,-2.061705) q[4];
u3(0.86548616,1.6366639,1.5738856) q[7];
cx q[1],q[7];
u3(1.0820356,-pi/2,pi/2) q[1];
u3(1.5750682,-3.0751615,3.0774699) q[7];
cx q[1],q[7];
u3(0.20615213,-pi,pi/2) q[1];
u3(1.5043758,0.0044345716,-1.6375118) q[7];
cx q[1],q[7];
u3(2.0754317,0.53341649,-1.9935856) q[1];
u3(2.7071316,0.42771754,1.2143663) q[7];
cx q[1],q[7];
u3(0.65895172,-pi/2,pi/2) q[1];
u3(1.5679319,-3.0750861,-3.0986131) q[7];
cx q[1],q[7];
u3(0.18414052,0,-pi/2) q[1];
u3(1.5043758,0.0044345716,-1.6375118) q[7];
cx q[1],q[7];
u3(1.3610174,-1.8691384,-0.93335) q[1];
u3(2.1173904,0.95262639,1.1804727) q[7];
u3(1.0952899,-2.2197553,-1.9542353) q[8];
cx q[6],q[8];
u3(0.54038152,-pi/2,pi/2) q[6];
u3(1.571682,-3.0750304,3.128307) q[8];
cx q[6],q[8];
u3(0.13589345,-pi,pi/2) q[6];
u3(1.5043758,0.0044345716,-1.6375118) q[8];
cx q[6],q[8];
u3(0.81931115,-1.8679372,-0.12680069) q[6];
cx q[4],q[6];
u3(0.67834443,-pi/2,pi/2) q[4];
u3(1.6099826,-3.0877667,2.5128906) q[6];
cx q[4],q[6];
u3(0.16274531,-pi,pi/2) q[4];
u3(1.5043758,0.0044345716,-1.6375118) q[6];
cx q[4],q[6];
u3(1.3863634,-2.5723099,-1.300467) q[4];
u3(1.1524053,2.0827076,1.2365809) q[6];
u3(1.8588745,-1.3886886,0.97465616) q[8];
u3(1.7660586,-2.6276978,0.67363081) q[9];
cx q[2],q[9];
u3(0.55998266,-pi/2,pi/2) q[2];
u3(1.5775761,-3.0753701,3.0397186) q[9];
cx q[2],q[9];
u3(0.23495822,0,-pi/2) q[2];
u3(1.5043758,0.0044345716,-1.6375118) q[9];
cx q[2],q[9];
u3(1.7107676,0.68012191,2.8466448) q[2];
cx q[2],q[8];
u3(0.84927149,-pi/2,pi/2) q[2];
u3(1.5887181,-3.0774789,2.8693876) q[8];
cx q[2],q[8];
u3(0.56505442,0,-pi/2) q[2];
u3(1.5043758,0.0044345716,-1.6375118) q[8];
cx q[2],q[8];
u3(1.4886029,1.0833189,-1.8751378) q[2];
cx q[3],q[2];
u3(1.6012858,-3.0824082,2.6664285) q[2];
u3(0.82749527,-pi/2,pi/2) q[3];
cx q[3],q[2];
u3(1.5043758,0.0044345716,-1.6375118) q[2];
u3(0.04173857,-pi,pi/2) q[3];
cx q[3],q[2];
u3(1.1079864,-0.44170064,-0.88908611) q[2];
u3(0.92699511,1.1815908,-0.52731818) q[3];
cx q[3],q[7];
u3(0.54647602,-pi/2,pi/2) q[3];
u3(1.5730105,-3.0750613,3.1083736) q[7];
cx q[3],q[7];
u3(0.072196264,-pi,pi/2) q[3];
u3(1.5043758,0.0044345716,-1.6375118) q[7];
cx q[3],q[7];
u3(1.3121869,0.31438522,-2.8988973) q[3];
u3(1.2493778,-2.4777294,-2.5560073) q[7];
cx q[3],q[7];
u3(0.94794036,-pi/2,pi/2) q[3];
u3(1.5821404,-3.0759968,2.9705914) q[7];
cx q[3],q[7];
u3(0.88943241,0,-pi/2) q[3];
u3(1.5043758,0.0044345716,-1.6375118) q[7];
cx q[3],q[7];
u3(1.4557097,-2.1164314,-2.9401935) q[3];
u3(1.9020028,2.2167133,-3.0127301) q[7];
u3(2.5944666,-0.61277619,2.5855595) q[8];
cx q[4],q[8];
u3(0.16079671,-pi/2,pi/2) q[4];
u3(1.5834452,-3.0762355,2.9506924) q[8];
cx q[4],q[8];
u3(0.061079909,-pi,pi/2) q[4];
u3(1.5043758,0.0044345716,-1.6375118) q[8];
cx q[4],q[8];
u3(0.34815574,-2.7610762,-2.099642) q[4];
u3(2.9077121,0.63226733,-2.3525407) q[8];
cx q[8],q[2];
u3(1.5770907,-3.0753223,3.0470353) q[2];
u3(1.2587789,-pi/2,pi/2) q[8];
cx q[8],q[2];
u3(1.5043758,0.0044345716,-1.6375118) q[2];
u3(0.15781923,-pi,pi/2) q[8];
cx q[8],q[2];
u3(1.4510567,0.47849017,-1.8957688) q[2];
u3(1.1336628,2.8434818,-0.42539904) q[8];
u3(1.9811377,-1.4316194,1.8082875) q[9];
cx q[0],q[9];
u3(1.0866218,-pi/2,pi/2) q[0];
u3(1.5738145,-3.0750928,3.0963051) q[9];
cx q[0],q[9];
u3(0.27047367,0,-pi/2) q[0];
u3(1.5043758,0.0044345716,-1.6375118) q[9];
cx q[0],q[9];
u3(0.21853601,0.35623631,-2.5299963) q[0];
cx q[6],q[0];
u3(1.5863879,-3.0768736,2.9055128) q[0];
u3(0.72830502,-pi/2,pi/2) q[6];
cx q[6],q[0];
u3(1.5043758,0.0044345716,-1.6375118) q[0];
u3(0.22709513,0,-pi/2) q[6];
cx q[6],q[0];
u3(2.6569258,1.4359672,1.0817192) q[0];
u3(1.5615501,2.6704889,2.1337115) q[6];
cx q[1],q[6];
u3(0.83903238,-pi/2,pi/2) q[1];
u3(1.569905,-3.0750304,-3.1282223) q[6];
cx q[1],q[6];
u3(0.016049311,-pi,pi/2) q[1];
u3(1.5043758,0.0044345716,-1.6375118) q[6];
cx q[1],q[6];
u3(0.78131334,2.978781,0.46744825) q[1];
u3(0.83720465,-0.0039112162,-3.0229704) q[6];
u3(2.4115639,-1.0292608,2.5636282) q[9];
cx q[5],q[9];
u3(1.1621769,-pi/2,pi/2) q[5];
u3(1.5781065,-3.0754265,3.0317166) q[9];
cx q[5],q[9];
u3(0.38793634,0,-pi/2) q[5];
u3(1.5043758,0.0044345716,-1.6375118) q[9];
cx q[5],q[9];
u3(0.27106001,1.5732428,2.5541892) q[5];
cx q[5],q[4];
u3(1.5832694,-3.0762018,2.9533775) q[4];
u3(0.58880305,-pi/2,pi/2) q[5];
cx q[5],q[4];
u3(1.5043758,0.0044345716,-1.6375118) q[4];
u3(0.024338894,-pi,pi/2) q[5];
cx q[5],q[4];
u3(1.8071395,2.5044034,0.38769861) q[4];
cx q[4],q[8];
u3(0.84400841,-pi/2,pi/2) q[4];
u3(2.0668185,2.0309484,-3.1278486) q[5];
cx q[5],q[1];
u3(1.5874121,-3.0771286,2.8896775) q[1];
u3(0.95163443,-pi/2,pi/2) q[5];
cx q[5],q[1];
u3(1.5043758,0.0044345716,-1.6375118) q[1];
u3(0.44158621,-pi,pi/2) q[5];
cx q[5],q[1];
u3(0.73811375,-0.93290923,1.311298) q[1];
cx q[1],q[3];
u3(0.82006025,-pi/2,pi/2) q[1];
u3(1.579471,-3.0755913,3.011099) q[3];
cx q[1],q[3];
u3(0.1522974,0,-pi/2) q[1];
u3(1.5043758,0.0044345716,-1.6375118) q[3];
cx q[1],q[3];
u3(1.0901878,-1.9151604,-2.8242536) q[1];
u3(1.8357801,-1.1578882,-3.0719351) q[3];
u3(1.6964946,-0.51460771,-0.93566608) q[5];
u3(1.5746391,-3.0751353,3.0839194) q[8];
cx q[4],q[8];
u3(0.010391532,-pi,pi/2) q[4];
u3(1.5043758,0.0044345716,-1.6375118) q[8];
cx q[4],q[8];
u3(1.480838,1.8456323,-3.1353026) q[4];
u3(2.3664543,1.7610336,0.90695396) q[8];
u3(1.6491748,2.7252539,2.453778) q[9];
cx q[0],q[9];
u3(0.87580537,-pi/2,pi/2) q[0];
u3(1.5705858,-3.0750248,-3.1384349) q[9];
cx q[0],q[9];
u3(0.3925893,-pi,pi/2) q[0];
u3(1.5043758,0.0044345716,-1.6375118) q[9];
cx q[0],q[9];
u3(2.4219245,-2.1335053,2.6144809) q[0];
cx q[0],q[2];
u3(0.55001101,-pi/2,pi/2) q[0];
u3(1.5787729,-3.0755034,3.0216539) q[2];
cx q[0],q[2];
u3(0.33858398,0,-pi/2) q[0];
u3(1.5043758,0.0044345716,-1.6375118) q[2];
cx q[0],q[2];
u3(0.7529846,1.8134073,-0.95537282) q[0];
cx q[0],q[8];
u3(0.93262376,-pi/2,pi/2) q[0];
u3(1.5413717,-1.8558806,-1.8025294) q[2];
cx q[5],q[2];
u3(1.5904154,-3.0779772,2.8428443) q[2];
u3(0.7803854,-pi/2,pi/2) q[5];
cx q[5],q[2];
u3(1.5043758,0.0044345716,-1.6375118) q[2];
u3(0.23126929,0,-pi/2) q[5];
cx q[5],q[2];
u3(2.5119878,-3.0334025,-2.0869122) q[2];
u3(0.41966851,2.3525694,-0.58009503) q[5];
u3(1.5697676,-3.0750324,-3.1261607) q[8];
cx q[0],q[8];
u3(0.21912678,0,-pi/2) q[0];
u3(1.5043758,0.0044345716,-1.6375118) q[8];
cx q[0],q[8];
u3(1.3411086,2.3594178,0.4685525) q[0];
u3(0.75012441,-0.66873406,-1.2661173) q[8];
cx q[8],q[3];
u3(1.5705232,-3.075025,-3.1374962) q[3];
u3(0.87333282,-pi/2,pi/2) q[8];
cx q[8],q[3];
u3(1.5043758,0.0044345716,-1.6375118) q[3];
u3(0.088881126,-pi,pi/2) q[8];
cx q[8],q[3];
u3(1.7087915,0.0025978891,0.26123665) q[3];
u3(1.3697264,0.3833777,-0.12636973) q[8];
u3(0.98712562,-2.7057049,1.6545382) q[9];
cx q[6],q[9];
u3(0.99969123,-pi/2,pi/2) q[6];
u3(1.5804436,-3.0757262,2.9963695) q[9];
cx q[6],q[9];
u3(0.061290784,-pi,pi/2) q[6];
u3(1.5043758,0.0044345716,-1.6375118) q[9];
cx q[6],q[9];
u3(2.9870974,-0.76507991,2.4319828) q[6];
cx q[6],q[4];
u3(1.5728399,-3.0750558,3.1109347) q[4];
u3(0.60011344,-pi/2,pi/2) q[6];
cx q[6],q[4];
u3(1.5043758,0.0044345716,-1.6375118) q[4];
u3(0.13173963,-pi,pi/2) q[6];
cx q[6],q[4];
u3(1.035239,2.2982357,0.21392332) q[4];
cx q[2],q[4];
u3(0.70076053,-pi/2,pi/2) q[2];
u3(1.566489,-3.0751638,-3.0769367) q[4];
cx q[2],q[4];
u3(0.0051549448,-pi,pi/2) q[2];
u3(1.5043758,0.0044345716,-1.6375118) q[4];
cx q[2],q[4];
u3(2.6652617,1.4163827,0.44222249) q[2];
u3(2.1248539,1.8141921,2.7630953) q[4];
cx q[2],q[4];
u3(0.63260606,-pi/2,pi/2) q[2];
u3(1.5682996,-3.0750713,-3.1041327) q[4];
cx q[2],q[4];
u3(0.20526619,0,-pi/2) q[2];
u3(1.5043758,0.0044345716,-1.6375118) q[4];
cx q[2],q[4];
u3(1.8993188,0.64241064,-0.32595269) q[2];
u3(2.2460818,1.0093182,0.70853764) q[4];
u3(2.5346615,-3.0846366,-1.6965167) q[6];
cx q[6],q[1];
u3(1.5783375,-3.0754524,3.0282301) q[1];
u3(0.71056458,-pi/2,pi/2) q[6];
cx q[6],q[1];
u3(1.5043758,0.0044345716,-1.6375118) q[1];
u3(0.16119549,-pi,pi/2) q[6];
cx q[6],q[1];
u3(1.6555887,-2.6736882,-0.41131002) q[1];
u3(2.0955483,-1.5443118,1.606551) q[6];
u3(2.5183511,-1.0278744,2.5125717) q[9];
cx q[7],q[9];
u3(1.3836276,-pi/2,pi/2) q[7];
u3(1.5735827,-3.0750827,3.0997841) q[9];
cx q[7],q[9];
u3(0.14929587,0,-pi/2) q[7];
u3(1.5043758,0.0044345716,-1.6375118) q[9];
cx q[7],q[9];
u3(2.247348,-1.1520568,-1.4939347) q[7];
cx q[7],q[0];
u3(1.5748793,-3.0751496,3.0803093) q[0];
u3(0.93195029,-pi/2,pi/2) q[7];
cx q[7],q[0];
u3(1.5043758,0.0044345716,-1.6375118) q[0];
u3(0.36803406,-pi,pi/2) q[7];
cx q[7],q[0];
u3(1.7081564,-1.117171,0.31110341) q[0];
cx q[0],q[3];
u3(0.99799358,-pi/2,pi/2) q[0];
u3(1.5847108,-3.0764929,2.9313155) q[3];
cx q[0],q[3];
u3(0.18868424,-pi,pi/2) q[0];
u3(1.5043758,0.0044345716,-1.6375118) q[3];
cx q[0],q[3];
u3(1.9539603,1.022228,1.3721691) q[0];
u3(2.1609722,2.2038965,-0.7801559) q[3];
cx q[0],q[3];
u3(1.093197,-pi/2,pi/2) q[0];
u3(1.5873954,-3.0771243,2.8899359) q[3];
cx q[0],q[3];
u3(0.68176503,0,-pi/2) q[0];
u3(1.5043758,0.0044345716,-1.6375118) q[3];
cx q[0],q[3];
u3(0.94133658,-1.4587803,-3.0721678) q[0];
u3(1.2093504,1.6478581,2.3601313) q[3];
u3(1.6898952,-1.3949247,-2.5955135) q[7];
cx q[7],q[8];
u3(1.0210636,-pi/2,pi/2) q[7];
u3(1.592349,-3.0786052,2.8123353) q[8];
cx q[7],q[8];
u3(0.59735408,0,-pi/2) q[7];
u3(1.5043758,0.0044345716,-1.6375118) q[8];
cx q[7],q[8];
u3(2.3632804,-2.7477163,-2.6679401) q[7];
u3(0.7572158,0.26015388,0.76160108) q[8];
u3(1.6271366,1.4945125,-3.1361557) q[9];
cx q[5],q[9];
u3(0.63015264,-pi/2,pi/2) q[5];
u3(1.5845451,-3.0764577,2.9338567) q[9];
cx q[5],q[9];
u3(0.035868723,-pi,pi/2) q[5];
u3(1.5043758,0.0044345716,-1.6375118) q[9];
cx q[5],q[9];
u3(2.0380566,-1.6465716,-0.93069253) q[5];
cx q[5],q[1];
u3(1.5972138,-3.0804837,2.7340505) q[1];
u3(0.56288939,-pi/2,pi/2) q[5];
cx q[5],q[1];
u3(1.5043758,0.0044345716,-1.6375118) q[1];
u3(0.29721964,0,-pi/2) q[5];
cx q[5],q[1];
u3(1.9065402,-0.42077756,-2.4190956) q[1];
cx q[1],q[7];
u3(0.91453965,-pi/2,pi/2) q[1];
u3(0.97636749,2.4061269,1.8344693) q[5];
cx q[5],q[2];
u3(1.5830778,-3.0761656,2.956303) q[2];
u3(1.1134156,-pi/2,pi/2) q[5];
cx q[5],q[2];
u3(1.5043758,0.0044345716,-1.6375118) q[2];
u3(0.024161938,0,-pi/2) q[5];
cx q[5],q[2];
u3(2.3377381,2.3871046,-0.37083252) q[2];
u3(0.83250271,2.2333952,-1.5931287) q[5];
u3(1.5748723,-3.0751492,3.0804148) q[7];
cx q[1],q[7];
u3(0.093148765,0,-pi/2) q[1];
u3(1.5043758,0.0044345716,-1.6375118) q[7];
cx q[1],q[7];
u3(1.4409432,-2.5557042,-1.2338122) q[1];
u3(2.2857222,2.9914021,-1.969548) q[7];
cx q[3],q[7];
u3(0.67249578,-pi/2,pi/2) q[3];
u3(1.5811461,-3.0758328,2.9857097) q[7];
cx q[3],q[7];
u3(0.1754893,0,-pi/2) q[3];
u3(1.5043758,0.0044345716,-1.6375118) q[7];
cx q[3],q[7];
u3(1.5543603,-2.9458229,0.34685069) q[3];
u3(1.1008797,-2.0499625,3.0271868) q[7];
u3(1.406456,1.1182803,-1.8283538) q[9];
cx q[9],q[6];
u3(1.5915042,-3.0783228,2.8257033) q[6];
u3(0.62531315,-pi/2,pi/2) q[9];
cx q[9],q[6];
u3(1.5043758,0.0044345716,-1.6375118) q[6];
u3(0.21284907,0,-pi/2) q[9];
cx q[9],q[6];
u3(1.09604,1.2339594,-1.4264221) q[6];
cx q[6],q[8];
u3(1.0238802,-pi/2,pi/2) q[6];
u3(1.5675535,-3.0751034,-3.092931) q[8];
cx q[6],q[8];
u3(0.33541983,-pi,pi/2) q[6];
u3(1.5043758,0.0044345716,-1.6375118) q[8];
cx q[6],q[8];
u3(2.027605,0.4905289,0.99730092) q[6];
cx q[5],q[6];
u3(0.73279643,-pi/2,pi/2) q[5];
u3(1.5808466,-3.0757864,2.9902565) q[6];
cx q[5],q[6];
u3(0.14951599,0,-pi/2) q[5];
u3(1.5043758,0.0044345716,-1.6375118) q[6];
cx q[5],q[6];
u3(0.52593515,-2.1919305,-0.43596344) q[5];
u3(2.4626635,0.2383953,2.3002696) q[6];
u3(1.7151432,-0.43803047,1.4339352) q[8];
cx q[8],q[1];
u3(1.5783375,-3.0754524,3.0282301) q[1];
u3(0.71056458,-pi/2,pi/2) q[8];
cx q[8],q[1];
u3(1.5043758,0.0044345716,-1.6375118) q[1];
u3(0.16119549,-pi,pi/2) q[8];
cx q[8],q[1];
u3(2.2657185,2.3009898,-2.1015574) q[1];
u3(0.85180418,-2.8162316,1.7169387) q[8];
u3(1.9326108,-1.378216,-1.0370707) q[9];
cx q[4],q[9];
u3(0.80731896,-pi/2,pi/2) q[4];
u3(1.5945146,-3.0793874,2.7777795) q[9];
cx q[4],q[9];
u3(0.5470182,0,-pi/2) q[4];
u3(1.5043758,0.0044345716,-1.6375118) q[9];
cx q[4],q[9];
u3(1.0798974,1.9604229,3.0892278) q[4];
cx q[4],q[2];
u3(1.5870468,-3.0770356,2.8953328) q[2];
u3(1.080097,-pi/2,pi/2) q[4];
cx q[4],q[2];
u3(1.5043758,0.0044345716,-1.6375118) q[2];
u3(0.24063227,0,-pi/2) q[4];
cx q[4],q[2];
u3(2.8708865,-1.2420325,3.1121085) q[2];
u3(1.4105338,-2.6884836,0.5501591) q[4];
u3(2.3961615,0.93188036,0.15500834) q[9];
cx q[9],q[0];
u3(1.5949745,-3.0795645,2.7703825) q[0];
u3(0.85859503,-pi/2,pi/2) q[9];
cx q[9],q[0];
u3(1.5043758,0.0044345716,-1.6375118) q[0];
u3(0.0590261,-pi,pi/2) q[9];
cx q[9],q[0];
u3(1.8504601,1.1299767,-0.084887301) q[0];
u3(2.7306591,-1.0678968,2.9771239) q[9];
Ideal simulator with sampling¶
!ouqu-tp ideal simulate --input-qasm-file=hoge.qasm --shots=100
0011100100
1001100100
0101111110
0000111000
1010000000
0011010110
0110111101
1111011010
1000010111
1101111101
0011000100
1001010010
0101111011
0010001011
0011111000
0110001101
0001110101
0111100100
0010001000
0101101010
0101010000
0000010111
0100111001
1010001101
1010001010
0010001111
1000100010
0111111011
1111110010
1111110101
0000110111
1000100011
0110000010
1111101101
1001110110
0111110010
1011101010
1000011001
1010111010
1010001011
0111110010
0000110011
0010100011
0010000110
0011000101
0011110011
1110011101
0111110111
1100001110
0110111101
0100001010
0000111101
0001001000
1011101010
1011101101
0000100001
1101011010
1001001010
1000000000
0010100110
0000110011
0011110101
1001110111
0110001110
0001010011
0011110101
0101110110
0110001111
0111001001
0000001011
0111100011
1011010100
1010000111
0011001110
1110111010
0100110001
1101011011
1100000100
1101000011
0110111101
1010010110
0100110001
1110111000
0011011100
0000000100
0110100000
1010000011
0000001000
1111110010
1000101011
1111001100
0101100100
0011000101
1001010110
1001110101
1110010110
0101110110
1010010100
0011100001
0010101000
Idela observable calculation¶
OpenFermion形式でobservableを記述する
observable_openfermion = """(-0.09706626816763086+0j) [] +
(-0.04530261550379928+0j) [X0 X1 Y2 Y3] +
(0.04530261550379928+0j) [X0 Y1 Y2 X3] +
(0.04530261550379928+0j) [Y0 X1 X2 Y3] +
(-0.04530261550379928+0j) [Y0 Y1 X2 X3] +
(0.1714128264477691+0j) [Z0] +
(0.16868898170361213+0j) [Z0 Z1] +
(0.1206252348339042+0j) [Z0 Z2] +
(0.16592785033770346+0j) [Z0 Z3] +
(0.1714128264477692+0j) [Z1] +
(0.16592785033770346+0j) [Z1 Z2] +
(0.1206252348339042+0j) [Z1 Z3] +
(-0.22343153690813464+0j) [Z2] +
(0.17441287612261575+0j) [Z2 Z3] +
(-0.22343153690813466+0j) [Z3]+
(-0.22343153690813466+0j) [Z9]"""
f = open('foo.of', 'w')
f.write(observable_openfermion)
f.close()
!ouqu-tp ideal getval --input-qasm-file=hoge.qasm --input-openfermion-file=foo.of
-0.1319437440019187
Estimation of observable with ideal sampling (with shot noise)¶
各パウリの項ごとに測定を行い期待値を評価する。
!ouqu-tp ideal sampleval --input-qasm-file=hoge.qasm --input-openfermion-file=foo.of --shots=100
-0.03719944803605122
!ouqu-tp ideal sampleval --input-qasm-file=hoge.qasm --input-openfermion-file=foo.of --shots=1000
-0.14253293055280988
!ouqu-tp ideal sampleval --input-qasm-file=hoge.qasm --input-openfermion-file=foo.of --shots=10000
-0.12848035163919114
!ouqu-tp ideal sampleval --input-qasm-file=hoge.qasm --input-openfermion-file=foo.of --shots=100000
-0.13094693306442934
!ouqu-tp ideal sampleval --input-qasm-file=hoge.qasm --input-openfermion-file=foo.of --shots=1000000
-0.13041773801162165
Noise simulator¶
State preparation, measurement, 1 qubit gate, 2 qubitゲートに対してノイズを導入して、上記のsimulationを実行。
prob = 0.01
shots = int(10)
pp= prob # preparation
pm = prob #measurement
p1 = prob # 1Q depolarizing noise
p2 = prob # 2Q depolarizing noise
Noisy sampling¶
!ouqu-tp noisy simulate --input-qasm-file=hoge.qasm --shots=$shots --p1=$p1 --p2=$p2 --pm=$pm --pp=$pp
1011111110
0110010010
1001001110
1100001000
1110100100
1111010011
1101001001
0010111111
0001010100
0101111010
Noisy observable calculation (without sampling)¶
!ouqu-tp noisy getval --input-qasm-file=hoge.qasm --input-openfermion-file=foo.of --p1=$p1 --p2=$p2 --pm=$pm --pp=$pp
-0.10143695200881075
Noisy sampling estimation¶
shots = 1000
!ouqu-tp noisy sampleval --input-qasm-file=hoge.qasm --input-openfermion-file=foo.of --shots=$shots --p1=$p1 --p2=$p2 --pm=$pm --pp=$pp
-0.11863618809924222
実機の制約条件を満たす量子回路へと変換¶
サンプルとして、グローバーと位相推定の両方を使う複雑な量子回路について考える。
# qubit_lists: QV回路を生成するための量子ビットサブセットのリストのリスト
nqubits =5
qubit_lists = [[i for i in range(nqubits)]]
# ntrials: サブセットごとに作成するランダム回路の数
ntrials = 1
import warnings
warnings.filterwarnings('ignore')
qv_circs, qv_circs_nomeas = qv.qv_circuits(qubit_lists, ntrials)
qv_circs_nomeas[0] = qiskit.compiler.transpile(qv_circs_nomeas[0], basis_gates=['u1','u2','u3','cx'])
qv_circs_nomeas[0][0].draw(output='mpl')
f = open('qv.qasm', 'w')
f.write(qv_circs_nomeas[0][0].qasm())
f.close()
実機制約を満たした回路へと変換
!ouqu-tp trance trance --input-qasm-file=qv.qasm --input-cnot-json-file=ouqu-tp/sample/created_Cnet.json >transpiled_qv.qasm
output_circuit =QuantumCircuit.from_qasm_file('transpiled_qv.qasm')
output_circuit.draw(output='mpl')
CNOTの代わりにCrossResonance gateに変換することも可能
!ouqu-tp trance trance_res --input-qasm-file=qv.qasm --input-cnot-json-file=ouqu-tp/sample/created_Cnet.json >transpiled_cr_qv.qasm
もう少し複雑な量子アルゴリズムのトランスパイル¶
import matplotlib.pyplot as plt
import numpy as np
import math
# importing Qiskit
import qiskit
from qiskit import QuantumCircuit, transpile, assemble, Aer
# import basic plot tools
from qiskit.visualization import plot_histogram
def example_grover_iteration():
"""Small circuit with 5/16 solutions"""
# 回路を作る
qc = QuantumCircuit(4)
# オラクル
qc.h([2,3])
qc.ccx(0,1,2)
qc.h(2)
qc.x(2)
qc.ccx(0,2,3)
qc.x(2)
qc.h(3)
qc.x([1,3])
qc.h(2)
qc.mct([0,1,3],2)
qc.x([1,3])
qc.h(2)
# 拡散
qc.h(range(3))
qc.x(range(3))
qc.z(3)
qc.mct([0,1,2],3)
qc.x(range(3))
qc.h(range(3))
qc.z(3)
return qc
# 制御されたグローバーの作成
grit = example_grover_iteration().to_gate()
grit.label = "Grover"
cgrit = grit.control()
def qft(n):
#n量子ビットのQFT回路を作成
"""Creates an n-qubit QFT circuit"""
circuit = QuantumCircuit(4)
def swap_registers(circuit, n):
for qubit in range(n//2):
circuit.swap(qubit, n-qubit-1)
return circuit
def qft_rotations(circuit, n):
#swapゲートなしで回路で最初にn量子ビットのQFTを実装する
"""Performs qft on the first n qubits in circuit (without swaps)"""
if n == 0:
return circuit
n -= 1
circuit.h(n)
for qubit in range(n):
circuit.cp(np.pi/2**(n-qubit), qubit, n)
qft_rotations(circuit, n)
qft_rotations(circuit, n)
swap_registers(circuit, n)
return circuit
qft_dagger = qft(4).to_gate().inverse()
qft_dagger.label = "QFT†"
# 量子回路の作成
t = 4 # カウントする量子ビットの数
n = 4 # 探索する量子ビットの数
qc = QuantumCircuit(n+t, t) # 古典ビットとn+t量子ビットの回路
# 全ての量子ビットを |+>に初期化する
for qubit in range(t+n):
qc.h(qubit)
# Begin controlled Grover iterations
iterations = 1
for qubit in range(t):
for i in range(iterations):
qc.append(cgrit, [qubit] + [*range(t, n+t)])
iterations *= 2
# Do inverse QFT on counting qubits
qc.append(qft_dagger, range(t))
# Measure counting qubits
qc.measure(range(t), range(t))
# Display the circuit
qc.draw()
┌───┐ »
q_0: ┤ H ├─────■────────────────────────────────────────────────────────────»
├───┤ │ »
q_1: ┤ H ├─────┼──────────■──────────■──────────────────────────────────────»
├───┤ │ │ │ »
q_2: ┤ H ├─────┼──────────┼──────────┼──────────■──────────■──────────■─────»
├───┤ │ │ │ │ │ │ »
q_3: ┤ H ├─────┼──────────┼──────────┼──────────┼──────────┼──────────┼─────»
├───┤┌────┴────┐┌────┴────┐┌────┴────┐┌────┴────┐┌────┴────┐┌────┴────┐»
q_4: ┤ H ├┤0 ├┤0 ├┤0 ├┤0 ├┤0 ├┤0 ├»
├───┤│ ││ ││ ││ ││ ││ │»
q_5: ┤ H ├┤1 ├┤1 ├┤1 ├┤1 ├┤1 ├┤1 ├»
├───┤│ Grover ││ Grover ││ Grover ││ Grover ││ Grover ││ Grover │»
q_6: ┤ H ├┤2 ├┤2 ├┤2 ├┤2 ├┤2 ├┤2 ├»
├───┤│ ││ ││ ││ ││ ││ │»
q_7: ┤ H ├┤3 ├┤3 ├┤3 ├┤3 ├┤3 ├┤3 ├»
└───┘└─────────┘└─────────┘└─────────┘└─────────┘└─────────┘└─────────┘»
c: 4/═══════════════════════════════════════════════════════════════════════»
»
« »
«q_0: ──────────────────────────────────────────────────────────────────»
« »
«q_1: ──────────────────────────────────────────────────────────────────»
« »
«q_2: ─────■────────────────────────────────────────────────────────────»
« │ »
«q_3: ─────┼──────────■──────────■──────────■──────────■──────────■─────»
« ┌────┴────┐┌────┴────┐┌────┴────┐┌────┴────┐┌────┴────┐┌────┴────┐»
«q_4: ┤0 ├┤0 ├┤0 ├┤0 ├┤0 ├┤0 ├»
« │ ││ ││ ││ ││ ││ │»
«q_5: ┤1 ├┤1 ├┤1 ├┤1 ├┤1 ├┤1 ├»
« │ Grover ││ Grover ││ Grover ││ Grover ││ Grover ││ Grover │»
«q_6: ┤2 ├┤2 ├┤2 ├┤2 ├┤2 ├┤2 ├»
« │ ││ ││ ││ ││ ││ │»
«q_7: ┤3 ├┤3 ├┤3 ├┤3 ├┤3 ├┤3 ├»
« └─────────┘└─────────┘└─────────┘└─────────┘└─────────┘└─────────┘»
«c: 4/══════════════════════════════════════════════════════════════════»
« »
« ┌───────┐┌─┐
«q_0: ─────────────────────────────────┤0 ├┤M├─────────
« │ │└╥┘┌─┐
«q_1: ─────────────────────────────────┤1 ├─╫─┤M├──────
« │ QFT† │ ║ └╥┘┌─┐
«q_2: ─────────────────────────────────┤2 ├─╫──╫─┤M├───
« │ │ ║ ║ └╥┘┌─┐
«q_3: ─────■──────────■──────────■─────┤3 ├─╫──╫──╫─┤M├
« ┌────┴────┐┌────┴────┐┌────┴────┐└───────┘ ║ ║ ║ └╥┘
«q_4: ┤0 ├┤0 ├┤0 ├──────────╫──╫──╫──╫─
« │ ││ ││ │ ║ ║ ║ ║
«q_5: ┤1 ├┤1 ├┤1 ├──────────╫──╫──╫──╫─
« │ Grover ││ Grover ││ Grover │ ║ ║ ║ ║
«q_6: ┤2 ├┤2 ├┤2 ├──────────╫──╫──╫──╫─
« │ ││ ││ │ ║ ║ ║ ║
«q_7: ┤3 ├┤3 ├┤3 ├──────────╫──╫──╫──╫─
« └─────────┘└─────────┘└─────────┘ ║ ║ ║ ║
«c: 4/═══════════════════════════════════════════╩══╩══╩══╩═
« 0 1 2 3 # 実行してみましょう。そして結果を確かめてみましょう。
aer_sim = Aer.get_backend('aer_simulator')
transpiled_qc = transpile(qc, aer_sim)
qobj = assemble(transpiled_qc)
job = aer_sim.run(qobj)
hist = job.result().get_counts()
plot_histogram(hist)
qiskit_quantum_counting = qiskit.compiler.transpile(qc, basis_gates=['u1','u2','u3','cx'])
f = open('quatum_counting.qasm', 'w')
f.write(qiskit_quantum_counting.qasm())
f.close()
input_circuit =QuantumCircuit.from_qasm_file('quatum_counting.qasm')
#input_circuit.draw()
!ouqu-tp trance trance --input-qasm-file=quatum_counting.qasm --input-cnot-json-file=ouqu-tp/sample/created_Cnet.json >transpiled_quantum_counting.qasm
output_circuit =QuantumCircuit.from_qasm_file('transpiled_quantum_counting.qasm')
output_circuit.draw()